Self-shielded gradient field coil for magnetic resonance imaging

ABSTRACT

A gradient coil for a magnetic resonance imaging apparatus ( 10 ) includes a primary coil ( 16 ) defining an inner cylindrical surface ( 60 ), and shield coil ( 18 ) or coils defining a coaxial outer cylindrical surface ( 62 ). Coil jumps ( 74 ) connect the primary and shield coils ( 16, 18 ). The coil jumps ( 74 ) define a non-planar current-sharing surface ( 64 ) extending between inner and outer contours ( 66, 68 ) that coincide with the inner and outer cylindrical surfaces ( 60, 62 ), respectively. The coil ( 16, 18, 74 ) defines a current path that passes across the current sharing surface ( 64 ) between the inner and outer contours ( 66, 68 ) a plurality of times. Optionally, some primary coil turns ( 70 ) are electrically interconnected to define an isolated primary sub coil (P 2 ) that together with a second shield (S 2 , S 2′ , S 2″ ) enables a discretely or continuously selectable field of view.

The following relates to the diagnostic imaging arts. It particularly relates to magnetic field gradient coils for use in short-bore magnetic resonance imaging scanners, and will be described with particular reference thereto. However, the following relates more generally to magnetic field gradient coils for various types of magnetic resonance imaging scanners, and to dedicated gradient coils for imaging reduced volumes of interest.

In a magnetic resonance imaging system, spatial encoding of the magnetic resonances is typically accomplished by generating magnetic field gradients in an imaging region. These gradients are produced by a set of magnetic field gradient coils. Gradient coils are typically provided for producing independent magnetic field gradients in each of the transverse x- and y-directions, and in the axial z-direction.

For various applications, such as producing gradients in short-bore magnetic resonance imaging scanners, or producing dedicated localized gradients over a reduced field of view (FoV) for imaging specific organs, it is desirable to employ relatively shorter gradient coils. In scaling gradient coils to reduced length while retaining high gradient peak amplitudes and linearity and bore diameter, coil turns become increasingly compressed near the coil extremities. Moreover, imaging distortions and artifacts due to eddy currents induced in a surrounding cold shield become more problematic. Furthermore, adaptive or variable field of view gradient coils are desired to adapt a given magnetic resonance scanner or magnetic field gradient coil set for a wider range of applications.

The present invention contemplates an improved apparatus and method that overcomes the aforementioned limitations and others.

According to one aspect, a gradient coil is disclosed for a magnetic resonance imaging apparatus. A primary coil defines an inner cylindrical surface. A shield coil defines an outer cylindrical surface coaxially aligned with the inner cylindrical surface and having a larger cylindrical radius than the inner cylindrical surface. A plurality of coil jumps electrically connect the primary and shield coils. The coil jumps define a non-planar current-sharing surface extending between an inner contour coinciding with the inner cylindrical surface and an outer contour coinciding with the outer cylindrical surface. The primary coil, shield coil , and coil jumps cooperatively define a current path that passes across the current-sharing surface between the inner and outer contours a plurality of times.

According to another aspect, a method is provided for producing a magnetic field gradient in a magnet bore of a magnetic resonance imaging apparatus. An electrical current is circulated through a primary coil that defines an inner cylindrical surface. The electrical current is circulated through a shield coil that defines an outer cylindrical surface coaxially aligned with the inner cylindrical surface and having a larger cylindrical radius than the inner cylindrical surface. The electrical current is communicated between the primary and shield coils via a plurality of non-planar coil jumps a plurality of times.

According to yet another aspect, a method is provided for designing a gradient coil for a magnetic resonance imaging apparatus. Current densities are computed on coaxial inner and outer cylindrical surfaces using constraints on the stored energy and gradient field linearity and uniformity. The computed current densities are generally non-zero at non-coplanar boundary contours that define coil edges. Primary coil turns, shield coil turns, and coil jumps are selected to approximate the computed current densities on the inner and outer cylindrical surfaces.

One advantage of designing the gradient coil with shared turns on the primary and shield coil is that it allows shorter shielded gradient coil designs.

Another advantage resides in reduced magnetic field gradient distortion due to eddy currents in the cold shield or other conductive elements of the magnetic resonance imaging scanner.

Another advantage resides in reduced current winding densities at the extremities of the gradient coil.

Yet another advantage resides in reduced I²R heating, particularly in regions of high winding density.

Numerous additional advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments.

The invention may take form in various components and arrangements of components, and in various process operations and arrangements of process operations. The drawings are only for the purpose of illustrating preferred embodiments and are not to be construed as limiting the invention.

FIG. 1 shows an exemplary magnetic resonance imaging apparatus, in partial cross-section, and auxiliary components for performing magnetic resonance imaging.

FIG. 2 shows an enlarged cross-sectional view of the magnetic resonance imaging apparatus which shows a flared arrangement of primary and shield gradient coils and a current-sharing surface for passing current back and forth between the primary and shield gradient coils.

FIG. 3 diagrammatically shows a sectional view of the geometrical arrangement of the primary and shield gradient coils of FIG. 2, with a suitable dimensional nomenclature superimposed.

FIG. 4 diagrammatically shows a perspective view of the geometrical arrangement of the primary and shield gradient coils.

FIG. 5 shows an unfolded half-quadrant layout of the X-gradient primary and shield gradient coils and the interconnecting coil jumps.

FIG. 6 shows exemplary isolated primary and shield coil turns of the X-gradient coil of FIG. 5.

FIG. 7 shows exemplary communicating primary and shield coil turns of the X-gradient coil of FIG. 5.

FIG. 8A shows a plot of gradient coil system energy versus the z-component of the current density at the edge of the primary coil for an exemplary transverse gradient coil design which has G_(x)=30 mT/m and the eddy current effect is—0.22%.

FIG. 8B shows a plot of gradient coil system energy versus residual eddy current effect for an exemplary transverse gradient coil design which has G_(x)=30 mT/m and J_(zp, edge)=17 kA/m.

FIG. 9 shows a plot of the z-component of the current density for the primary and shield coils layout shown in FIG. 5.

FIG. 10 shows a plot of X-gradient non-uniformity at a radius of 0.001 meters for the primary and shield coils layout shown in FIG. 5.

FIG. 11 shows a plot of X-gradient non-linearity at z=0 meters as a function of radial position for the primary and shield coils layout shown in FIG. 5.

FIG. 12 shows a plot of residual eddy current distribution as a function of axial position Z produced by the coil layout shown in FIG. 5.

FIG. 13 shows a plot comparing eddy current effect as a function of axial position Z for the continuous current of FIG. 9 and the discretized coils of FIG. 5.

FIG. 14 shows an unfolded half-quadrant layout of a coil P2 formed by selecting fourteen turns on the primary coil of FIG. 5 that do have jumps.

FIG. 15 shows an unfolded half-quadrant layout of a coil P1 formed by disconnecting fourteen isolated primary coil turns of the primary gradient coil P2 of FIG. 14.

FIG. 16 shows an unfolded half-quadrant layout of a second shield coil S2 that operates in conjunction with the coil P1 to define a second field of view.

FIG. 17 shows a plot of X-gradient non-uniformity for the two discrete fields of view of the gradient coil illustrated in FIGS. 5 and 14-16.

FIG. 18 shows a plot of X-gradient non-linearity for the two discrete fields of view of the gradient coil illustrated in FIGS. 5 and 14-16.

FIG. 19 shows an unfolded half-quadrant layout of a second shield coil S2′ that operates in conjunction with the coils of FIGS. 5, 14, and 15 to define a discretely variable gradient coil with three fields of view.

FIG. 20 shows two gradient coils driven in parallel, with nomenclature specified for computing a slew rate for the parallel-combined coils.

FIG. 21 shows a suitable configuration of two primary coils located on a single surface and two shield coils designed as described herein that provide a continuously variable field of view.

FIG. 22 shows a plot of the z-component of the current density for primary and shield coils providing a 32 cm DSV field of view.

FIG. 23 shows a plot of the z-component of the current density for primary and shield coils of a correction or shim coil that cooperates with the 32 cm DSV coil of FIG. 22 to produce a continuously variable field of view.

FIG. 24 shows an unfolded half-quadrant layout of the primary and shield gradient coils and the interconnecting coil jumps of the discretized current densities of the 32 cm DSV coil of FIG. 22.

FIG. 25 shows an unfolded half-quadrant layout of the primary gradient coil of the correction coil.

FIG. 26 shows an unfolded half-quadrant layout of the shield gradient coil of the correction coil.

With reference to FIG. 1, a magnetic resonance imaging scanner 10 includes a cylindrical main magnet 12, which is preferably superconducting and cryoshielded. The main magnet 12 defines a magnet bore 14 inside of which a patient or other imaging subject is placed for imaging. The main magnet 12 produces a spatially and temporally constant and uniform main magnetic field oriented along a longitudinal axis of the bore 14. Instead of a superconducting magnet, a non-superconducting magnet can be used. Moreover, a vertical magnet, an open magnet, or other type of main magnet can be employed instead of the illustrated horizontal cylindrical main magnet 12. In one preferred embodiment, the main magnet 12 is a short-bore magnet.

Magnetic field gradient coils include a primary gradient coil 16 and a shield gradient coil 18 that cooperatively produce magnetic field gradients in the bore 14 for spatially encoding magnetic resonance signals, for producing magnetization-spoiling field gradients, or the like. Preferably, the magnetic field gradient coils 16, 18 include coils configured to produce magnetic field gradients in three orthogonal directions including transverse x- and y-directions. In addition to the shield coil 18, a cold shield 20 provides a high conductivity eddy current surface for residual gradient fields thus protecting the magnet coils still further away. In a suitable embodiment, the cold shield 20 is integrated into a housing of the main magnet 12.

A whole body radio frequency coil assembly 22 generates radio frequency pulses for exciting magnetic resonances. The radio frequency coil assembly 22 also serves to detect magnetic resonance signals. Optionally, additional local radio frequency coils or phased radio frequency coil arrays (not shown) are included for exciting and/or detecting magnetic resonances at localized areas in the bore 14.

Gradient pulse amplifiers 30 deliver controlled electrical currents to the magnetic field gradient coils 16, 18 to produce selected magnetic field gradients. A radio frequency transmitter 34, preferably digital, applies radio frequency pulses or pulse packets to the radio frequency coil assembly 22 to generate selected magnetic resonance excitations. A radio frequency receiver 36 also coupled to the radio frequency coil assembly 22 receives magnetic resonance signals. If more than one radio frequency coil is provided (such as a local coil or phased coil array), then different coils are optionally used for the magnetic resonance excitation and detection operations.

To acquire magnetic resonance imaging data of a subject, the subject is placed inside the magnet bore 14, preferably at or near an isocenter of the main magnetic field. A sequence controller 40 communicates with the gradient amplifiers 30 and the radio frequency transmitter 34 to produce selected transient or steady state magnetic resonance configurations in the subject, to spatially encode such magnetic resonances, to selectively spoil magnetic resonances, or otherwise generate selected magnetic resonance signals characteristic of the subject. The generated magnetic resonance signals are detected by the radio frequency receiver 36, and stored in a k-space memory 42. The imaging data is reconstructed by a reconstruction processor 44 to produce an image representation that is stored in an image memory 46. In one suitable embodiment the reconstruction processor 44 performs an inverse Fourier transform reconstruction.

The resultant image representation is processed by a video processor 48 and displayed on a user interface 50, which is preferably a personal computer, workstation, or other type of computer. Rather than producing a video image, the image representation can be processed by a printer driver and printed, transmitted over a computer network or the Internet, or the like. Preferably, the user interface 50 also allows a radiologist or other operator to communicate with the magnetic resonance sequence controller 40 to select magnetic resonance imaging sequences, modify imaging sequences, execute imaging sequences, and so forth.

With reference to FIG. 2-4, the magnet 12 of the magnetic resonance imaging scanner 10 is preferably a short-bore magnet defining a cylindrical magnet bore 14. The primary and shield gradient coils 16, 18 coincide with and define axially concentric cylindrical surfaces 60, 62, respectively. The inner and outer concentric cylindrical surfaces 60, 62 have radii R_(p), R_(s), respectively. The shield gradient coil 18 extends axially a half-length L_(s) along the outer cylindrical surface 62, while the primary gradient coil 16 extends axially a half-length L_(p) along the inner cylindrical surface 60. The cold shield 20 has a length 2 L_(CS). Preferably, the gradient coils 16, 18 and the cold shield 20 are bilaterally symmetric along the axial direction, and L_(s), L_(p), L_(CS) designate half-lengths.

To provide reduced length, increased openness, and improved shielding, the axial half-length L_(s) of the shield 30 gradient coil 18 is preferably greater than the axial half-length L_(p) of the primary gradient coil 16. The two gradient coils 16, 18 are connected along a non-planar current-sharing surface 64 that extends between an inner contour 66 coinciding with the inner cylindrical surface 60, and an outer contour 68 coinciding with the outer cylindrical surface 62. In a preferred embodiment, the non-planar current-sharing surface corresponds to a curved surface of a frustum of a cone, where the cone has an angle a respective to an axis of the axially concentric cylindrical surfaces 60, 62. The angle a corresponds to an arctangent of a ratio of a difference between the cylindrical surface radii R_(p), R_(s) and a difference between the primary and secondary or shield gradient coil half-lengths L_(p), L_(s). However, the non-planar current-sharing surface can also have a parabolic, piece-wise conic, or other shape.

With continuing reference to FIG. 2-4 and with further reference to FIG. 5, the primary gradient coil 16 includes a fingerprint-like pattern of conductive windings or current paths coinciding with the inner cylindrical surface 60, while the shield gradient coil 18 includes a fingerprint-like pattern of conductive windings or current paths coinciding with the outer cylindrical surface 62.

With continuing reference to FIGS. 2-5 and with further reference to FIGS. 6 and 7, the primary gradient coil 16 includes isolated primary coil turns 70 that connect with other primary gradient coil turns on the surface 60. The primary gradient coil 16 also includes communicating primary coil turns 72 that each connect with at least one coil jump 74 that coincides with the current-sharing surface 64. Similarly, the shield gradient coil 18 includes isolated shield coil turns 80 that connect with other shield coil turns on the surface 62, as well as communicating shield coil turns 84 that each connect with at least one coil jump 74.

As shown in FIG. 6, each isolated primary or shield coil turn 70, 80 defines a current loop on the cylindrical surface 60, 62, respectively. As shown in FIG. 7, the communicating primary and shield coil turns 72, 82 and two coil jumps 74 define a current loop that passes back and forth across the current-sharing surface 64 via the two coil jumps 74.

Those skilled in the art will recognize that the gradient coils 16, 18 have substantial advantages over past coils, especially for short-magnet designs. By connecting the primary coil 16 and secondary or shield coil 18 with a plurality of coil jumps 74, high coil turn densities near the coil edges are reduced. Moreover, by axially flaring the shield coil 18 relative to the primary coil 16 improved shielding is obtained versus a design using primary and shield coils of equal lengths. The axial flaring of the shield coil 18 is accommodated by the non-planar current-sharing surface 64 to permit a combination of a flared shield and multiple coil jumps.

With returning reference to FIG. 3, a suitable method for designing the gradient coils 16, 18 uses the following parameters. The angle α is called a flare angle. The axial coil half-lengths L_(p), L_(s), and the cold shield half-length L_(CS) are shown and were discussed previously. Preferably, L_(p)<L_(s)<L_(CS), with L_(CS) substantially larger than L_(s). The primary coil 16 coincides with and defines the inner cylindrical surface 60 which has a radius R_(p). The shield coil 18 coincides with and defines the outer cylindrical surface 62 which has a radius R_(s). The cold shield 20 is preferably cylindrical with a radius R_(CS), and is preferably coaxial with the cylindrical surfaces 60, 62. For a preferred embodiment, R_(p)<R_(s)<R_(CS).

The gradient coil design goals include producing: (i) given gradient strength G_(X) or G_(Y) at the isocenter of the magnet bore 14; (ii) high slew rate; (iii) good quality gradient magnetic field characteristics, such as linearity and uniformity, within a given diameter spherical volume (DSV); and (iv) good shielding. It is recognized that leakage of the magnetic field into the conductive cold shield 20 will be present. This leakage is responsible for the residual eddy currents induced in the cold shield 20 when the gradient coil 16, 18 is operated. These residual eddy currents can affect the quality of the gradient field inside the DSV and are characterized by the residual eddy current effect, which is defined as follows: $\begin{matrix} {\left. {{ed\_ ef}(r)} \right|_{DSV} = \left. \frac{B_{z}^{eddy\_ current}(r)}{B_{z}^{Grad}(r)} \middle| {}_{DSV}{*100\quad\%} \right.} & (1) \end{matrix}$ where the symbol |_(DSV) indicates that the residual eddy current effect is evaluated at the surface of the DSV. One skilled in is the art knows that if the above characteristic is reasonably constant as a function of spatial position, r, then the residual eddy current effect can be substantially corrected with gradient pre-emphasis. Thus, the absolute level of eddy current magnitude, although generally desired or preferred to be as small as possible, is not necessarily the most important aspect.

In designing a self-shielded gradient coil, which from the beginning has finite lengths of the primary and shield coils, shielding conditions for the normal component of the total magnetic field B_(r)(r) on the surface of the shield coil are typically imposed according to: $\begin{matrix} \begin{matrix} \begin{matrix} {\left. {B_{\rho}(r)} \right|_{\rho = R_{x}} = 0} \\ {{{ɛ(k)} = \left. {{F_{s}(k)} + {{R(k)}{F_{p}(k)}}}\rightarrow 0 \right.},\quad{{R(k)} = \frac{R_{p}{I_{1}^{\prime}\left( {kR}_{p} \right)}}{R_{s}{I_{1}^{\prime}\left( {kR}_{s} \right)}}}} \end{matrix} \\ {{F_{p,s}(k)} = {\overset{L_{p,s}}{\int\limits_{- L_{p,s}}}{{f_{\varphi}^{({p,s})}(z)}{\cos({kz})}{\mathbb{d}z}}}} \end{matrix} & (2) \end{matrix}$ where f_(φ) ^((p,s))(z) is the φ-component of the current density on the primary (p) or shield (s) coil, R_(p,s), L_(p,s) are the radii and half-length of the primary (p) or shield (s) coil, respectively, and I₁′(kR) is a derivative of a Bessel function of the first kind I₁(kR). Equation (2) holds for any k, and is substantially restrictive. The shielding error function ε(k) is a measure of how well the gradient coil is shielded. To design a well-shielded gradient coil Equation (2) leads to high current density on the primary coil with a negative lobe near its edge, leading therefore to high inductance.

If the residual eddy current effect is non-zero it should have small variation over the imaging volume to be correctable. This means, that the variation of the residual eddy current effect should be small, even if the residual eddy current effect is about, say, 10%. Therefore, the coil is designed to satisfy the conditions that there is no high current density near the ends of the coils and the residual eddy current effect is correctable.

The gradient coil design includes the following process operations: (i) allowing the current to flow off the end of the primary coil 16 to the shield coil 18 along the conical or other non-planar current-sharing surface 64 disposed therebetween; (ii) the cold shield 20 is included in the coil design; and (iii) the residual eddy current effect is considered as a variable in the design together with gradient field quality characteristics. The latter will be introduced as point constraints.

In the following, an exemplary X-gradient coil design is described. Those skilled in the art can readily adapt the exemplary X-gradient coil design to the designing of a Y-gradient coil. The current distribution on the X-gradient coil in terms of cylindrical coordinates has the following form: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{J(r)} = {{e_{\rho}{J_{\rho}(r)}} + {e_{\varphi}{J_{\varphi}(r)}} + {e_{z}{J_{z}(r)}}}} \\ {{J_{\rho}(r)} = {{J^{(\rho)}(\rho)}\left\{ {{\delta\left( {z - {Z(\rho)}} \right)} + {\delta\left( {z + {Z(\rho)}} \right)}} \right\}{\Theta\left( {\rho - R_{p}} \right)}{\Theta\left( {R_{s} - \rho} \right)}{\sin(\varphi)}}} \end{matrix} \\ {{J_{\varphi}(r)} = \left\{ {{{f_{\varphi}^{(p)}(z)}{\delta\left( {\rho - R_{p}} \right)}{\Theta\left( {z + L_{p}} \right)}{\Theta\left( {L_{p} - z} \right)}} +} \right.} \end{matrix} \\ {\left. {{f_{\varphi}^{(s)}(z)}{\delta\left( {\rho - R_{s}} \right)}{\Theta\left( {z + L_{s}} \right)}{\Theta\left( {L_{s} - z} \right)}} \right\}{\cos(\varphi)}} \end{matrix} \\ {{J_{z}(r)} = \left\{ {{{f_{z}^{(p)}(z)}{\delta\left( {\rho - R_{p}} \right)}{\Theta\left( {z + L_{p}} \right)}{\Theta\left( {L_{p} - z} \right)}} + {{f_{z}^{(s)}(z)}{\delta\left( {\rho - R_{s}} \right)}\Theta}} \right.} \\ {{\left. {\left( {z + L_{s}} \right){\Theta\left( {L_{s} - z} \right)}} \right\}{\sin(\varphi)}} + {{J^{(z)}(\rho)}\left\{ {{\delta\left( {z - {Z(\rho)}} \right)} - {\delta\left( {z + {Z(\rho)}} \right)}} \right\}\Theta}} \\ {\left( {\rho - R_{p}} \right){\Theta\left( {R_{s} - \rho} \right)}{\sin(\varphi)}} \\ {{{Z(\rho)} = {L_{p} + {\left( {\rho - R_{p}} \right){\cot(\vartheta)}}}},\quad{{\cot(\vartheta)} = \frac{L_{s} - L_{p}}{R_{s} - R_{p}}}} \end{matrix} & (3) \end{matrix}$ The continuity equation ∇·J(r)=0 applied to the current J(r) given in Equation (3) above leads to the following equations: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{{R_{p}{f_{z}^{{(p)}^{\prime}}(z)}} = {f_{\varphi}^{(p)}(z)}},\quad{{R_{s}{f_{z}^{{(s)}^{\prime}}(z)}} = {f_{\varphi}^{(s)}(z)}},} \\ {{\frac{\partial\left( {\rho\quad{J^{(\rho)}(\rho)}} \right)}{\partial\rho} = {\left. 0\rightarrow{J^{(\rho)}(\rho)} \right. = \frac{C}{\rho}}},} \end{matrix} \\ {{{J^{(\rho)}\left( R_{p} \right)} = {f_{z}^{(p)}\left( L_{p} \right)}},{{J^{(\rho)}\left( R_{s} \right)} = {- {f_{z}^{(s)}\left( L_{s} \right)}}},} \end{matrix} \\ {{{J^{(z)}(\rho)} = {{J^{(\rho)}(\rho)}{\cot(\vartheta)}}},} \end{matrix} \\ {C = {{R_{p}{f_{z}^{(p)}\left( L_{p} \right)}} = {{- R_{s}}{f_{z}^{(s)}\left( L_{p} \right)}}}} \end{matrix} & (4) \end{matrix}$ where the functions f_(z) ^(p,s)(z) are suitably represented as Fourier series according to: $\begin{matrix} {{{f_{z}^{(p)}(z)} = {\sum\limits_{n = 1}^{\quad}{a_{n}^{(p)}{\sin\left( {\frac{\pi\quad n}{2L_{p}}z} \right)}}}},{{f_{z}^{(s)}(z)} = {\sum\limits_{n = 1}^{\quad}{a_{n}^{(s)}{\sin\left( {\frac{\pi\quad n}{2L_{s}}z} \right)}}}}} & (5) \end{matrix}$ with the coefficients α_(n) ^((p)), α_(n) ^((s)) determined through extremization of the following functional: $\begin{matrix} {W = {E - {\sum\limits_{i}^{\quad}{\left( {{B_{z}^{coil}\left( r_{i} \right)} - R_{i}} \right)\lambda_{i}}} - {\sum\limits_{j}^{\quad}{\left( {{B_{z}^{CS}\left( r_{j} \right)} - {C_{j}{B_{z}^{coil}\left( r_{j} \right)}}} \right)\Lambda_{j}}} - {\left( {{f_{z}^{p}\left( L_{p} \right)} - A_{p}} \right)\Lambda_{p}} - {\left( {{f_{z}^{s}\left( L_{s} \right)} + {\frac{R_{p}}{R_{s}}A_{p}}} \right)\Lambda_{s}}}} & (6) \end{matrix}$ where: E represents total energy stored in the magnetic fields of the gradient coil 16, 18 and the cold shield 20; B_(z) ^(coil)(r_(i)) is the z-component of the coil magnetic field at a set of points r_(i) inside the DSV where the field has predetermined values B_(i); B_(z) ^(CS)(r_(j)) is the z-component of the magnetic field produced by the cold shield 20 at a set of points r_(j) on the surface of the DSV; B_(z) ^(coil)(r_(j)) is the z-component of the coil's magnetic field at the same set of points; and C_(j) is a set of numbers that govern the residual eddy current effect. Finally, λ_(i) and Λ_(j) are the corresponding Lagrange multipliers. The quantity A_(p) is a desired value of the z-component of the current density on the primary coil at its edge, and the quantity R_(p)/R_(s)A_(p) is the value of the z-component of the current density at the edge of the shield coil. The continuity conditions are represented by the last two terms in Equation (6) and Λ_(p), Λ_(s) are the corresponding Lagrange multipliers for the current continuity. The quantities A_(p) and C_(j) are not variables and can be chosen arbitrarily. In deriving the Equation (6), we have assumed that the normal component of the total magnetic field vanishes on the surface of the magnet cold shield 20, which has radius R_(CS), which is assumed much longer than the length of the gradient set. This gives the following equation for the cold shield residual eddy current: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{f_{\varphi}^{CS}(z)} = {\frac{1}{\pi}{\int\limits_{0}^{\infty}{{F_{\varphi}^{CS}(k)}{\cos({kz})}{\mathbb{d}k}}}}} \\ {{F_{\varphi}^{CS}(k)} = {{{- {R\left( {k,R_{p},R_{CS}} \right)}}{F_{\varphi}^{p}(k)}} - {{R\left( {k,R_{s},R_{CS}} \right)}{F_{\varphi}^{s}(k)}} -}} \end{matrix} \\ {{2\frac{R_{p}{f_{z}^{p}\left( L_{p} \right)}}{R_{CS}}{\overset{R_{s}}{\int\limits_{R_{p}}}{\frac{I_{1}\left( {k\quad\rho} \right)}{I_{1}^{\prime}\left( {kR}_{CS} \right)}{\cos\left( {{kZ}(\rho)} \right)}\frac{\mathbb{d}\rho}{\rho}}}},} \end{matrix} \\ {{R\left( {k,R,R_{CS}} \right)} = \frac{{RI}_{1}^{\prime}({kR})}{R_{CS}{I_{1}^{\prime}\left( {kR}_{CS} \right)}}} \end{matrix} & (7) \end{matrix}$ where the Fourier transforms F_(p,s)(k) of the primary and the shield current distributions are given in Equation (2). Variation of the functional W of Equation (6) with respect to α_(n) ^(p), α_(n) ^(s) determines the solution in terms of continuous current distributions.

Optimization of the functional given in Equation (6) with respect to α_(n) ^((p)) and α_(n) ^((s)) determines a suitable solution by adjusting the value of the current density at the edge of the primary coil. After a solution is found the continuous current distributions on primary coil, shield coils and the interconnecting surface are discretized such that each loop on the primary coil carries the same current I, and each loop on the shield coil carries opposite current. The discretized coil is tested for its performance: gradient strength at isocenter, linearity and uniformity of the gradient strength over the DSV, and slew rate and residual eddy current effect.

The continuous current distributions are discretized and approximated by coils 70, 72, 74, 80, 82 carrying current along a current path. A suitable fingerprint pattern of the primary and shield coils 16, 18 and coil jumps 74 determined by the forgoing process is shown in FIG. 5. The left and right half of the coil are the primary and shield coils 16, 18 respectively joined by coil jumps 74 at fixed angular positions φ along the current-sharing surface 64 that connects the concentric inner and outer cylindrical surfaces 60, 62.

With continuing reference to FIG. 5 and with further reference to FIGS. 6 and 7, the two kinds of current paths, namely the isolated primary and shield coil turns 70, 80 exemplarily shown in FIG. 6 and the communicating primary and shield coil turn 72, 82 shown in FIG. 7 are generated. The isolated current loops 70, 80 reside only on the primary or shield coil 16, 18, that is, only on one of the two cylindrical surfaces 60, 62. The communicating current paths 72, 82 define current loops that reside on both the primary and shield coils 16, 18, that is, on both of the two cylindrical surfaces 60, 62, with coil jumps 74 defining current jumps or interconnections between them.

The B_(z) field from all of these paths is calculated to determine residual eddy currents in the cold shield 20, which can be compared to the continuous current residual eddy current effects on the DSV. Each isolated loop 70, 80 resides on either the primary coil 16 or shield coil 18 and is closed on itself. Examples of these loops are shown in FIG. 6, which shows an isolated primary loop 70 on the inner cylindrical surface 60 and an isolated shield coil loop 80 located in the first quadrant, for example. Each loop carries a current I. The current on the isolated primary coil loop 70 flows clockwise, while the current on the shield coil loop 80 flows counterclockwise. Exemplary communicating primary and shield loops 72, 82 defining a current loop passing back and forth across the non-planar current sharing surface 64 is shown in FIG. 7.

A suitable process for calculating the magnetic field and associated residual eddy current in the magnet cold shield 20 from the discretized coil 16, 18 is the path-integral approach. The coil is discretized, so an equation describing each loop is known, i.e. the equations φ=f_(p)(z), φ=f_(s)(z). For the isolated loops shown in FIG. 6, the conditions f_(p)(p1)=f_(p)(p2)=0, f_(s)(s1)=f(s2)=0 apply. For the communicating loops shown in FIG. 7, different conditions f_(p)(p1)=f_(s)(s1)=f_(p)f(L_(p))=φ_(f), f_(s)(L_(s))=φ_(f) apply, where φ_(f) is the angle at which the current-sharing coil jump occurs. Rather than the parameters p₁, p₂, s₁, s₂, the description henceforth uses z_(i) and z_(f).

The current distribution produced by the communicating loops 72, 82 shown in FIG. 7 is given by: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{J(r)} = {{e_{\rho}{J_{\rho}(r)}} + {e_{\varphi}{J_{\varphi}(r)}} + {e_{z}{J_{z}(r)}}}} \\ {{{J_{\rho}(r)} = {{J^{\rho}(\rho)}{\delta\left( {z - {Z(\rho)}} \right)}{\Theta\left( {\rho - R_{p}} \right)}{\Theta\left( {R_{s} - \rho} \right)}\left\{ {{\delta\left( {\varphi - \varphi_{f}} \right)} - {\delta\left( {\varphi + \varphi_{f}} \right)}} \right\}}},} \end{matrix} \\ {{J_{\varphi}(r)} = {Q_{p}R_{p}{f_{p}^{\prime}(z)}\left\{ {{\delta\left( {\varphi - {f_{p}(z)}} \right)} + {\delta\left( {\varphi + {f_{p}(z)}} \right)}} \right\}{\delta\left( {\rho - R_{p}} \right)}\Theta}} \end{matrix} \\ {{\left( {z_{f}^{p} - z} \right){\Theta\left( {z - z_{i}^{z}} \right)}} + {Q_{s}R_{s}{f_{s}^{\prime}(z)}\left\{ {{\delta\left( {\varphi - {f_{s}(z)}} \right)} + {\delta\left( {\varphi + {f_{s}(z)}} \right)}} \right\}{\delta\left( {\rho - R_{s}} \right)}\Theta}} \end{matrix} \\ {{\left( {z_{f}^{s} - z} \right){\Theta\left( {z - z_{i}^{s}} \right)}},} \end{matrix} \\ {{J_{z}(r)} = {{B_{p}\left\{ {{\delta\left( {\varphi - {f_{p}(z)}} \right)} - {\delta\left( {\varphi + {f_{p}(z)}} \right)}} \right\}{\delta\left( {\rho - R_{p}} \right)}{\Theta\left( {z_{f}^{p} - z} \right)}{\Theta\left( {z - z_{i}^{p}} \right)}} +}} \end{matrix} \\ {{B_{s}\left\{ {{\delta\left( {\varphi - {f_{s}(z)}} \right)} - {\delta\left( {\varphi + {f_{s}(z)}} \right)}} \right\}{\delta\left( {\rho - R_{s}} \right)}{\Theta\left( {z_{f}^{s} - z} \right)}{\Theta\left( {z - z_{i}^{s}} \right)}} +} \end{matrix} \\ {{J^{z}(\rho)}{\delta\left( {z - {Z(\rho)}} \right)}{\Theta\left( {\rho - R_{p}} \right)}{\Theta\left( {R_{s} - \rho} \right)}\left\{ {{\delta\left( {\varphi - \varphi_{f}} \right)} - {\delta\left( {\varphi + \varphi_{f}} \right)}} \right\}} \end{matrix} & (8) \end{matrix}$ where in Equation (8) Q_(p), B_(p), Q_(s), B_(s) are unknown constants and J^(ρ)(ρ), J^(z)(ρ) are unknown functions to be determined from the continuity equation ∇·J(r)=0. Thus, $\begin{matrix} \begin{matrix} {{Q_{p} = {B_{p} = \frac{I}{R_{p}}}},} \\ {{Q_{s} = {B_{s} = {- \frac{I}{R_{s}}}}},} \\ {{{J^{z}(\rho)} = {{J^{\rho}(\rho)}{\cot(\vartheta)}}},} \\ {{J^{\rho}(\rho)} = {\frac{I}{\rho}.}} \end{matrix} & (9) \end{matrix}$ This allows expressions for all components of the vector potentials, and therefore for magnetic field and residual eddy current, to be found. For example, the expression for the z-component of the magnetic field produced by a loop, combined from all four quadrants, is given by: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{{B_{z}(r)} = {{B_{z}^{p}(r)} + {B_{z}^{s}(r)} + {B_{z}^{CJ}(r)}}},} \\ {{B_{s}^{p}(r)} = {{- \frac{8\mu_{o}I}{\pi^{2}}}{\sum\limits_{n = 1}^{\infty}{{\cos\left( {\left( {{2n} - 1} \right)\varphi} \right)}{\int\limits_{0}^{\infty}{{\mathbb{d}k}\quad{\cos({kz})}\left( {kR}_{p} \right){{\mathfrak{J}}_{{2n} - 1}^{p}(k)}}}}}}} \end{matrix} \\ {\left\{ {{\Theta\left( {R_{p} - \rho} \right)}{I_{{2n} - 1}\left( {k\quad\rho} \right)}{K_{{2n} - 1}^{\prime}\left( {kR}_{p} \right)}{\Theta\left( {\rho - R_{p}} \right)}{K_{{2n} - 1}\left( {k\quad\rho} \right)}{I_{{2n} - 1}^{\prime}\left( {kR}_{p} \right)}} \right\},} \end{matrix} \\ {{B_{z}^{s}(r)} = {\frac{8\mu_{o}I}{\pi^{2}}{\sum\limits_{n = 1}^{\infty}{{\cos\left( {\left( {{2n} - 1} \right)\varphi} \right)}{\int\limits_{0}^{\infty}{{\mathbb{d}k}\quad\cos\quad({kz})\left( {kR}_{s} \right){{\mathfrak{J}}_{{2n} - 1}^{s}(k)}}}}}}} \end{matrix} \\ {\left\{ {{\Theta\left( {R_{s} - \rho} \right)}{I_{{2n} - 1}\left( {k\quad\rho} \right)}{K_{{2n} - 1}^{\prime}\left( {kR}_{s} \right)}{\Theta\left( {\rho - R_{s}} \right)}{K_{{2n} - 1}\left( {k\quad\rho} \right)}{I_{{2n} - 1}^{\prime}\left( {kR}_{s} \right)}} \right\},} \end{matrix} \\ {{B_{z}^{CJ}(r)} = {{- \frac{8\mu_{o}I}{\pi^{2}}}{\sum\limits_{n = 1}^{\infty}{\left( {{2n} - 1} \right){\cos\left( {\left( {{2n} - 1} \right)\varphi} \right)}{\sin\left( {\left( {{2n} - 1} \right)\varphi_{f}} \right)}}}}} \end{matrix} \\ {\int\limits_{0}^{\infty}{{\mathbb{d}k}\quad{\cos({kz})}{\int\limits_{R_{p}}^{R_{s}}{\frac{\mathbb{d}x}{x}{\cos\left( {{kZ}(x)} \right)}\left\{ {{{\Theta\left( {\rho - x} \right)}{I_{{2n} - 1}\left( {k\quad\rho} \right)}{K_{{2n} - 1}({kx})}} +} \right.}}}} \end{matrix} \\ \left. {{\Theta\left( {X - \rho} \right)}{I_{{2n} - 1}({kx})}{K_{{2n} - 1}\left( {k\quad\rho} \right)}} \right\} \end{matrix} \\ {{where}\text{:}} \end{matrix} & (10) \\ \begin{matrix} {{{{\mathfrak{J}}_{{2n} - 1}^{p}(k)} = {\overset{z_{f}^{p}}{\int\limits_{z_{i}^{p}}}{{\mathbb{d}{{zf}_{p}^{\prime}(z)}}{\cos\left\lbrack {\left( {{2n} - 1} \right){f_{p}(z)}} \right\rbrack}{\cos({kz})}}}},} \\ {{{\mathfrak{J}}_{{2n} - 1}^{s}(k)} = {\overset{z_{f}^{s}}{\int\limits_{z_{i}^{s}}}{{\mathbb{d}{{zf}_{s}^{\prime}(z)}}{\cos\left\lbrack {\left( {{2n} - 1} \right){f_{s}(z)}} \right\rbrack}{{\cos({kz})}.}}}} \end{matrix} & (11) \end{matrix}$ It will be recognized that all odd modes in cos((2n−1)φ) are present. These expressions hold for both isolated and communicating loops. For isolated coils 70, 80 that are closed on themselves, z_(i) ^(p,s), z_(f) ^(p,s) are the initial and final points of the primary (p) or shield (s) loop in the upper-half of the first quadrant. Communicating loops 72, 82 which connect with at least one coil jump 74 have z_(i) ^(p,s) as initial points on the primary (p) or shield (s) loop in the upper half of the first quadrant, and z_(f) ^(p,s)=L_(p,s). From Equation (2) the result for the z-component of the residual eddy current is: $\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{{f_{z}^{CS}\left( {z,\varphi} \right)} = {\sum\limits_{n = 1}^{\infty}{{f_{z,{{2n} - 1}}^{CS}(z)}{\cos\left\lbrack {\left( {{2n} - 1} \right)\varphi} \right\rbrack}}}},} \\ {{{f_{z,{{2n} - 1}}^{CS}(z)} = {\frac{1}{\pi}{\int\limits_{0}^{\infty}{{F_{{2n} - 1}^{CS}(k)}{\sin({kz})}{\mathbb{d}z}}}}},} \end{matrix} \\ {{F_{{2n} - 1}^{CS}(k)} = {{- \frac{8I}{\pi}}\left( {{\frac{1}{R_{p}}{\Re\left( {k,R_{p},R_{CS}} \right)}{\aleph_{{2n} - 1}^{p}(k)}} - {\frac{1}{R_{s}}{\Re\left( {k,R_{s},R_{CS}} \right)}}} \right.}} \end{matrix} \\ {{\aleph_{{2n} - 1}^{s}(k)} + {\frac{\sin\left\lbrack {\left( {{2n} - 1} \right)\varphi_{f}} \right.}{{kR}_{CS}}\left\lbrack {{\frac{R_{CS}}{R_{p}}{\Re\left( {k,R_{p},R_{CS}} \right)}{\cos\left( {kz}_{f}^{p} \right)}} -} \right.}} \end{matrix} \\ {{\frac{R_{CS}}{R_{s}}{\Re\left( {k,R_{s},R_{CS}} \right)}{\cos\left( {kz}_{f}^{s} \right)}} + {\left( {{2n} - 1} \right)^{2}{\overset{R_{s}}{\int\limits_{R_{p}}}{\frac{\mathbb{d}x}{x}{\cos\left( {{kZ}(x)} \right)}}}}} \end{matrix} \\ \left. \left. \frac{I_{{2n} - 1}({kx})}{\left( {kR}_{CS} \right){I_{{2n} - 1}^{\prime}\left( {kR}_{CS} \right)}} \right\rbrack \right) \end{matrix} & (12) \end{matrix}$ where the following notations are introduced: $\begin{matrix} \begin{matrix} {{{\Re\left( {k,R,R_{CS}} \right)} = {\left( \frac{R_{p}}{R_{s}} \right)^{2}\frac{I_{{2n} - 1}^{\prime}({kR})}{I_{{2n} - 1}^{\prime}\left( {kR}_{CS} \right)}}},} \\ {{\aleph_{{2n} - 1}^{p,s}(k)} = {\overset{z_{f}^{p,s}}{\int\limits_{z_{i}^{p,s}}}{{\mathbb{d}z}\quad{\sin\left\lbrack {\left( {{2n} - 1} \right){f_{p,s}(z)}} \right\rbrack}{{\sin({kz})}.}}}} \end{matrix} & (13) \end{matrix}$

The above approach can also be adapted to the design of a self-shielded axial gradient (Z-gradient) coil. In this case, it is sufficient to find the J_(φ)component of the current density on both the primary and the shield coils. This component does not necessarily vanish at the coil edges. In other words, the current continuity equation ∇·J=0 does not imply that the J_(φ) component of the current density should be a continuous function of z even at the edges of the coils. Therefore, the values of the current densities at the edges of the coils can be treated as free parameters during the design. Thus, J(r)=e _(φ) J _(φ)(r) J _(φ)(r)={f _(φ) ^((p))(z)δ(ρ−R _(p))Θ)(z+L _(p))Θ(L _(p) −z)+f _(φ) ^((s))(z)δ(ρ−R_(s))Θ(z+L _(s))Θ(L _(s) −z)}  (14) where the functions f_(φ) ^(p,s)(z) are suitably represented as a Fourier series: $\begin{matrix} {{{f_{\varphi}^{(p)}(z)} = {\sum\limits_{n = 1}^{\quad}{a_{n}^{(p)}{\sin\left( {\frac{\pi\quad n}{2L_{p}}z} \right)}}}},{{f_{\varphi}^{(s)}(z)} = {\sum\limits_{n = 1}^{\quad}{a_{n}^{(s)}{\sin\left( {\frac{\pi\quad n}{2L_{s}}z} \right)}}}}} & (15) \end{matrix}$ with the coefficients α_(n) ^((p)), α_(n) ^((s)) being determined from extremization of the functional: $\begin{matrix} {W = {E - {\sum\limits_{i}^{\quad}{\left( {{B_{z}^{coil}\left( r_{i} \right)} - B_{i}} \right)\lambda_{i}}} - {\sum\limits_{j}^{\quad}{\left( {{B_{z}^{CS}\left( r_{j} \right)} - {C_{j}{B_{z}^{coil}\left( r_{j} \right)}}} \right)\Lambda_{j}}} - {\left( {{f_{z}^{p}\left( L_{p} \right)} - A_{p}} \right)\Lambda_{p}} - {\left( {{f_{z}^{s}\left( L_{s} \right)} + A_{s}} \right)\Lambda_{s}}}} & (16) \end{matrix}$ where: E is the total energy of the gradient coil and the cold shield; B_(z) ^(coil)(r_(i)) is the z-component of the coil magnetic field at a set of points r_(i) inside the DSV where the field has predetermined values B_(i); B_(z) ^(CS)(r_(j)) is the z-component of the magnetic field produced by the cold shield eddy currents at a set of points r_(j) on the surface of the DSV; B_(z) ^(coil)(r_(j)) is the z-component of the coil's magnetic field at the same set of points; and C_(j) is a set of numbers that govern the residual eddy current effect. Finally, λ_(i) and Λ_(j) are the corresponding Lagrange multipliers. The quantities A_(p), A_(s) are the desired values of the φ-component of the current densities at the edges of the primary and shield coils, respectively. In Equation (16) the shielding condition on the surface of the cold shield 20, namely: B _(ρ)(r)|ρ=R _(CS) =0   (17) has been incorporated. This leads to a continuous residual eddy current distribution given by: $\begin{matrix} \begin{matrix} {{f_{\varphi}^{CS}(z)} = {\frac{1}{\pi}{\int\limits_{0}^{\infty}{{F_{CS}(k)}{\sin({kz})}{\mathbb{d}z}}}}} \\ {{F_{CS}(k)} = {{{- {{Sh}\left( {k,R_{p},R_{CS}} \right)}}{F_{p}(k)}} - {{{Sh}\left( {k,R_{s},R_{CS}} \right)}{F_{s}(k)}}}} \\ {{{F_{p,s}(k)} = {\overset{L_{p,s}}{\int\limits_{- L_{ps}}}{{f_{\varphi}^{p,s}(z)}{\sin({kz})}{\mathbb{d}z}}}},} \\ {{{Sh}\left( {k,R,R_{CS}} \right)} = {\frac{{RI}_{1}({kR})}{R_{CS}{I_{1}\left( {kR}_{CS} \right)}}.}} \end{matrix} & (18) \end{matrix}$ Optimization of the functional given in Equation (16) with respect to α_(n) ^((p)) and α_(n) ^((s)) determines a suitable solution by adjusting the values of the current densities at the edge of the primary coil and the shield coil. After a solution is found the continuous current distributions on both primary and shield coils are discretized such that each loop on the primary coil carries the same current I, and each loop on the shield coil carries opposite current. The discretized coil is tested for its performance: gradient strength at isocenter, linearity and uniformity of the gradient strength over the DSV, slew rate, and residual eddy current effect.

An exemplary transverse gradient coil design is described. The exemplary coil is designed for a magnet bore having a length of 1.2 meters and a bore radius of 0.45 meters. Design dimensions of the exemplary transverse gradient coil are as follows: R_(P)=0.34385 m, L_(P)=0.434 meters, R_(S)=0.435224 meters, L_(S)=0.555 meters. In all examples described herein the target eddy current effect was chosen to be −0.22%.

The coil is aligned such that the flare angle α shown in FIG. 3 is 36.9°. In designing the X-gradient coil the field constraints shown in Table 1 were used. TABLE 1 ρ-position Φ- z-position (m) position (m) B_(z) (Tesla) 0.001 0.000 0.000 2.7e−5 0.180 0.000 0.000 4.86528530650126e−3 0.001 0.000 0.200 2.069587125e−5

The first constraint corresponds to the gradient strength of G_(x)=27 mT/m at the isocenter, the second constraint is assigned for no-rollover of the z-component of the gradient magnetic field, and the third constraint corresponds to −23% of non-uniformity of the gradient strength at z=0.2 meters. The residual eddy current effect constraints were imposed at a set of points, which are given in Table 2. TABLE 2 z-position (m) ρ-position (m) 0.0128 0.2495 0.06304 0.2373 0.1296 0.1904 0.19904 0.0245 The radius of the cold shield 20 was chosen to be R_(CS)=0.475 m. The coil jumps 74 between primary and shield coils 16, 18 at the ends of the gradient assembly result in a significant reduction in stored energy and a change of residual eddy current effects. A study was done of total energy vs. the current density at the edge of the primary coil (quantity A_(p) in Equation 6) and the residual eddy current effect (quantity C_(j), the percent measurement of the residual eddy current effect, in Equation 6).

FIG. 8A shows the plot for the energy of the system, rescaled to the gradient strength of 30 mT/m and corrected back to account for the loss of the gradient strength due to the residual eddy current effect at the isocenter of the coil, as a function of the current density at the edge of the primary coil when the residual eddy current effect is equal to an exemplary −0.22%. FIG. 8B, shows the plot for the energy of the system, rescaled to the gradient strength of 30 mT/m and corrected back to account for the loss of the gradient strength due to the residual eddy current effect at the isocenter of the coil, as a function of the eddy current effect. In FIG. 8B the current density at the edge of the primary coil is equal to an exemplary 17.0 kA/m. FIG. 8A shows that the total energy at the gradient strength of G_(x)=30 mT/m after accounting for the residual eddy current effect at the isocenter of the coil, has a minimum at a certain value of the quantity A_(p) for a given value of the residual eddy current effect. This value is suitably adjusted such that when the current densities on the primary and shield coils 16, 18 are discretized, the numbers of turns are integers.

FIG. 9 shows the results for the continuous current distributions shown in the unfolded coil layout of FIG. 5. The continuous current distributions shown in FIG. 9 were discretized with 29 turns on the primary coil 16 and 19 turns on the shield coil 18 per quadrant with 14 coil jumps 74 per half of each quadrant. Each closed loop carries current I=398.64 A in order to produce gradient strength G_(x)=30.0 mT/m at isocenter. The current layout per half quadrant is shown in FIG. 5. Henceforth, P and S denote the whole primary coil 16 and the whole shield coil 18, respectively, which together with the coil jumps 74 provide 50 cm×40 cm imaging volume.

The discretized current was used to calculate magnetic field quality characteristics. FIG. 10 shows X-gradient field non-uniformity (taking the transverse gradient coil as producing an X-gradient; of course, the transverse gradient coil can also be suitably oriented to produce a Y-gradient or other transverse gradient), while FIG. 11 shows X-gradient field non-linearity.

FIG. 12 shows the z-component of residual eddy current produced by the X-gradient coil having the current paths shown in FIG. 5. The residual eddy current of FIG. 12 and the discretized loops shown in FIG. 5 were used to calculate a residual eddy current effect over 50 cm×40 cm volume. The comparison of this result with that obtained from the continuous currents is shown in FIG. 13. This figure shows that the discretization procedure does not introduce substantial error in terms of magnitude and its variation over the volume.

With reference to FIGS. 14 and 15, the transverse gradient coil whose layout is shown in FIG. 5 is readily adapted to provide a discretely variable field of view. In a suitable modification, fourteen isolated primary coils 70 near the eye of the primary coil fingerprint are electrically decoupled from the remainder to define an interconnected coil shown in FIG. 14 and designated as P2 herein. As shown in FIG. 15, a remainder P1 of the primary coil 16 remains electrically connected with the shield coil 18. For convenience, the original coupled primary coil 16 is designated as P herein. Hence, P=P1+P2. Similarly, the shield coil 18 is also designated by S herein.

With reference to FIG. 16, a second shield coil S2 is designed to compensate for the fourteen electrically disconnected loops comprising the primary coil P2. A radial position of the second shield coil S2 is 0.415224 m. The current layout on the second shield coil is shown in FIG. 16. The current on the second shield S2 is I=398.64 A. Combination of the coil shown in FIG. 15 with the second shield coil shown in FIG. 16, where the current flows in the same direction as on the primary coil, provides a reduced DSV of 30 cm. The non-uniformity and non-linearity of the gradient strength of two modes gradient coils are shown in FIGS. 17 and 18, respectively. The main mode coil corresponding to the layout of FIG. 5 provides 50 cm×40 cm imaging volume, while the coil of FIG. 15 combined with the second shield S2 shown in FIG. 16 provides 30 cm DSV.

If a second power supply is provided, the gradient coil 16, 18 can be further modified to provide three discrete fields of view. The second shield coil S2 is replaced by an eight-turn shield coil S2′ positioned at the same radial location. A base mode, which provides 50 cm×40 cm imaging volume, is generated by the gradient coil shown in FIG. 5 driven at a current I.

A second mode, which provides 30 cm DSV, is generated by the gradient coil shown in FIG. 15 driven with a current I, combined with the second shield S2 shown in FIG. 19 driven with a current I/2. Hence, two power supplies (or a current divider) are used.

A third mode is generated by the gradient coil shown in FIG. 15 operating concurrently with the inner fourteen loop coil P2 of FIG. 14. The current in the primary coil P2 is reversed relative to the current in the primary coil P1 of FIG. 15. The second shield S2 of FIG. 19 is driven with a current −I. In this third mode, all coils are driven in series, and one power supply is used. This configuration provides 26 cm DSV.

When two coils are driven in series, a slew rate of the series-combined coils is given by: $\begin{matrix} {{SR} = {S*\frac{V - {IR}}{L}}} & (19) \end{matrix}$ where: S is a sensitivity of the coil: S=G/I, with G and I being the gradient strength and the current, respectively; V is a output voltage of the power supply; R is a resistance of the coil including the cable resistance; and L is an inductance of the coil including the cable and filters inductance.

With reference to FIG. 20, when two coils are driven in parallel, the slew rate of the combined coil can be found by considering the coupling equations between two electric circuits, given by: $\begin{matrix} {{{{L_{1}\frac{\mathbb{d}I_{1}}{\mathbb{d}t}} + {R_{1}I_{1}} + {M\frac{\mathbb{d}I_{2}}{\mathbb{d}t}}} = V_{1}},{{{L_{2}\frac{\mathbb{d}I_{2}}{\mathbb{d}t}} + {R_{2}I_{2}} + {M\frac{\mathbb{d}I_{1}}{\mathbb{d}t}}} = V_{2}}} & (20) \end{matrix}$ where L₁, R₁, V₁ and L₂, R₂, V₂ are the inductance, resistance, and output voltage of the first and the second circuits, respectively, and M is the mutual inductance between these two coils. Then the slew rate of the combined coil is given by: $\begin{matrix} \begin{matrix} {{SR} = {\frac{1}{\Delta}\left\{ {{\left( {{S_{1}L_{2}} - {S_{2}M}} \right)\left( {V_{1} - {R_{1}I_{1}}} \right)} + {\left( {{S_{2}L_{1}} - {S_{1}M}} \right)\left( {V_{2} - {R_{2}I_{2}}} \right)}} \right\}}} \\ {\Delta = {{L_{1}L_{2}} - M^{2}}} \end{matrix} & (21) \end{matrix}$ where S₁, S₂ is the sensitivity of the first and second coils, respectively.

Moreover, the gradient coils of FIG. 5, including the primary coils 16 and secondary or shield coils 18, which are also designated herein as P and S, respectively, are readily modified to provide a continuous field of view.

Starting with the initial coil system P+S, which provides the imaging volume with largest FoV, a variable FoV system can be achieved in the following way: (1) electrically disconnect several turns near the eye of the primary coil (this defines the coil P2, for example as shown in FIG. 14) and the remaining primary coil is P1=P−P2 (for example as shown in FIG. 15). The coil P1 includes communicating coil turns 72 that communicate with communicating coil turns 82 of the shield coil 18 by coil jumps 74 at the ends of the gradient. The positions of these interconnects are known and fixed. Optionally, the coil P1 also has some isolated turns that are only on the inner cylindrical surface 60.

A second shield S1 is designed to conform with the following conditions: (i) the number of conical jumps on this coil and their angular and axial positions match those on the primary coil P1; and (ii) shield S1 together with the primary coil P1 form a well-shielded pair. The system P1+S1 provides a reduced DSV1<DSV relative to the gradient coil 16, 18 of FIG. 5.

A second shield coil S2″ is designed for the coil P2. The second shield coil S2″ together with the coil P2 form a well shielded pair. The combination of S1 and S2″ is equivalent to the initial coil S. The original coil P is not used (it was an intermediate process operation in the design procedure). Rather, two well-shielded pairs of coils: P1+S1 and P2+S2″ are the final result of the design.

These coil sets can be used as to provide either a discretely selectable FoV as in the previous examples or a continuously variable FoV. For a discretely selectable FoV, the coil (P1+P2)+(S1+S2″) provides the initial imaging volume DSV. The coil P1+S1 provides DSV1<DSV. The coil (P1−P2)+(S1−S2″) provides DSV2<DSV1.

With reference to FIG. 21, for a continuously variable FoV, the coil sets (P1+S1) and (P2+S2″) provide a continuously variable FoV when used in conjunction with two power supplies. In this arrangement, one coil set (P1+S1) is driven by a current I from a first power supply, while the coil set (P2+S2″) is driven by the current a*I from a second power supply. The parameter α adjusts the size of the FoV.

With reference to FIG. 22, another gradient coil with a continuously variable field of view is described. In this approach, a well-shielded coil is designed in accordance with process operations described previously, including current transfer between inner and outer cylindrical surfaces corresponding to primary and shield gradient coils, respectively, across a non-planar current-sharing surface. FIG. 22 shows a suitable continuous current density for such a coil which produces a 32 cm DSV.

A correction or shim gradient coil is designed in such a way that it produces zero gradient strength at an isocenter of the coil and together with the 32 cm DSV gradient coil produces a 50 cm×40 cm DSV. By itself the correction coil is well shielded.

The correction coil is suitably designed as follows: a well-shielded gradient coil is designed, including a non-planar current-sharing surface, that provides 50 cm×40 cm imaging volume and has the same gradient strength at the isocenter as that of the 32 cm DSV coil. The value of the current density at the edge of the primary coil is the same as in the gradient coil that provides 32 cm DSV. The difference (subtraction) between the current densities (on primaries and separately shield coils) of 50 cm×40 cm and 32 cm DSV defines the correction coil. The values of the current density at the edge of the coils are adjusted in such a way that allows proper discretization. An example of the current density on the correction gradient coil is shown in FIG. 23.

The current distributions of the 32 cm DSV gradient coil are discretized to generate a current layout such as that shown in FIG. 24. Similarly, the correction or shim coil current distributions are discretized to generate a current layout such as the primary and shield shim coil current layouts shown in FIGS. 25 and 26, respectively. Variation of the amount of current in the correction gradient coil allows continuous variation of the FoV. The driving coil interconnections of FIG. 21 can be used in driving the 32 cm FOV gradient coils of FIG. 24 and the correction coils of FIGS. 25 and 26. In this approach the coils P1 and S1 correspond to the well-shielded coil set of FIG. 24 that provides the 32 cm FoV, and P2 and S2″ correspond to the well-shielded correction coil pair shown in FIGS. 25 and 26.

FIGS. 24-26 describe an arrangement in which the primary coil (FIG. 24) produces a base field of view (32 cm) that is continuously enlargeable by cooperative operation of the correction coils of FIGS. 25 and 26. Alternatively, a base coil can provide a large base field of view which is continuously reducible by cooperative operation of correction coils.

Moreover, although design of transverse X-gradient coils are described herein, those skilled in the art recognize that the design of a transverse Y-gradient coil is substantively similar to design of an X-gradient coil, differing by a rotation of the whole gradient assembly by 90°.

The gradient coil design procedures described herein are readily modified to accommodate other design aspects. For example, local forces on the gradient coil can be reduced by introducing an additional constraint that the net (thrust) force on the gradient coil be zero. This constraint incorporates information about the main magnet that the gradient coil will be used with. Reduction in the local forces on the gradient coil can reduce gradient tube deflection when the coil is operated.

The invention has been described with reference to the preferred embodiments. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof. 

1. A gradient coil for a magnetic resonance imaging apparatus, the gradient coil including: a primary coil defining an inner cylindrical surface; a shield coil defining an outer cylindrical surface coaxially aligned with the inner cylindrical surface and having a larger cylindrical radius than the inner cylindrical surface; and a plurality of coil jumps electrically connecting the primary and shield coils, the coil jumps defining a non-planar current-sharing surface extending between an inner contour coinciding with the inner cylindrical surface and an outer contour coinciding with the outer cylindrical surface; the primary coil, shield coil, and coil jumps cooperatively defining a current path that passes across the current-sharing surface between the inner and outer contours a plurality of times.
 2. The gradient coil as set forth in claim 1, wherein: the primary coil extends axially a primary coil length along the inner cylindrical surface; and the shield coil extends axially a shield coil length along the outer cylindrical surface; the shield coil length not equal to the primary coil length.
 3. The gradient coil as set forth in claim 2, wherein: the current sharing surface corresponds to a curved surface of a frustum of a cone with a cone angle defined by a difference between radii of the inner and outer cylindrical surfaces and a difference between the primary and shield coil lengths.
 4. The gradient coil as set forth in claim 1, wherein: the inner and outer contours are circular contours, and the current-sharing surface includes a conical surface portion extending between the inner and outer circular contours.
 5. The gradient coil as set forth in claim 1, wherein the primary coil includes: communicating primary coil turns that electrically connect with a coil jump; and isolated primary coil turns that do not electrically connect with a coil jump.
 6. The gradient coil as set forth in claim 5, wherein the shield coil includes: communicating shield coil turns that electrically communicate with communicating primary coil turns via connecting coil jumps.
 7. The gradient coil as set forth in claim 5, wherein at least some of the isolated primary coil turns are electrically interconnected to define an isolated primary sub-coil, and the gradient coil further includes: a switch having at least: a first state in which the isolated primary sub-coil is electrically connected with the communicating primary coil turns, and a second state in which the isolated primary sub-coil is electrically isolated from the communicating primary coil turns; the first and second states corresponding to first and second selectable fields of view.
 8. The gradient coil as set forth in claim 7, wherein the isolated primary sub-coil is deenergized in the second state.
 9. The gradient coil as set forth in claim 7, wherein the isolated primary sub-coil is energized with opposite polarities in the two states.
 10. The gradient coil as set forth in claim 7, wherein the gradient coil further includes: a second shield coil that is energized in one of the two states to improve uniformity of the corresponding field of view.
 11. The gradient coil as set forth in claim 5, wherein at least some isolated primary coil turns are interconnected to define a selectively electrically switched primary sub-coil, the gradient coil further including: a second shield coil that is selectively energized in conjunction with switching of the primary sub-coil to define a variable field of view.
 12. The gradient coil as set forth in claim 1, further including: a shielded correction coil that cooperatively adjusts a field of view over a continuous range.
 13. The gradient coil as set forth in claim 1, further including: a generally cylindrical cold shield coaxially aligned with the outer cylindrical surface and having a larger cylindrical radius than the outer cylindrical surface, the cold shield carrying eddy current that produces a substantially spatially constant residual eddy current effect.
 14. The gradient coil as set forth in claim 13, wherein the substantially spatially constant residual eddy current effect is non-zero.
 15. The gradient coil as set forth in claim 1, wherein the gradient coil is a transverse gradient coil.
 16. The gradient coil as set forth in claim 1, wherein the coil jumps are selected to minimize the stored energy of the coil.
 17. A magnetic resonance scanner comprising: a main magnet for generating a temporally constant magnetic field; a gradient coil as set forth in claim 1 for inducing magnetic field gradients across the temporally constant magnetic field; at least one RF coil disposed adjacent the gradient coil; an RF transmitter connected with one of the RF coils for inducing and manipulating resonance; an RF receiver connected with one of the RF coils for demodulating induced resonance; and a reconstruction processor for reconstructing the demodulated resonance into an image representation.
 18. A method for producing a magnetic field gradient in a magnet bore of a magnetic resonance imaging apparatus, the method including: circulating an electrical current through a primary coil that defines an inner cylindrical surface; circulating the electrical current through a shield coil that defines an outer cylindrical surface coaxially aligned with the inner cylindrical surface and having a larger cylindrical radius than the inner cylindrical surface; and communicating the electrical current back and forth between the primary and shield coils via a plurality of non-planar coil jumps a plurality of times.
 19. The method as set forth in claim 18, further including: selecting a total number of coil jumps which minimizes the stored energy of the gradient coil.
 20. The method as set forth in claim 19, wherein the selecting of a total number of coil jumps includes: constraining selection of the coil jumps to produce substantially spatially constant residual eddy current effect from a cold shield that surrounds the shield coil.
 21. The method as set forth in claim 18, further including: computing current densities on the inner and outer cylindrical surfaces using constraints including minimizing stored energy and minimizing the variation of the residual eddy current effect, the current densities being generally non-zero at the inner and outer contours; arranging coil turns of the primary and shield coils to approximate the computed current densities on the inner and outer cylindrical surfaces; and during the arranging of coil turns arranging coil jumps to approximate the computed non-zero current densities at the inner and outer contours.
 22. The method as set forth in claim 21, further including: simultaneously with the computing of current densities on the inner and outer cylindrical surfaces, computing current densities on a current-sharing surface, the arranging of coil jumps being further constrained to approximate the computed current densities on the current-sharing surface.
 23. The method as set forth in claim 21, wherein the computing of current densities further includes: constraining the current densities to produce a substantially spatially constant eddy current effect produced by a current density in a cold shield that surrounds the shield coil.
 24. The method as set forth in claim 23, further including: applying a gradient pre-emphasis that substantially cancels a magnetic field produced by the eddy currents in the cold shield.
 25. The method as set forth in claim 18, wherein the primary coil (16) includes a plurality of coil loops, the method further including: selectively electrically isolating at least some primary coil loops from the communicating of the electrical current between the primary and shield coils via the plurality of coil jumps, the selective electrical isolating defining a coil set which combined with a second shield coil provides a second selectable field of view.
 26. The method as set forth in claim 18, wherein the primary coil includes a plurality of coil loops, the method further including: selectively removing at least some primary coil loops from the communicating of the electrical current between the primary and shield coils via the plurality of coil jumps, the selectively removed primary coil loops being interconnected to define a primary sub-coil; and circulating second electrical current through the primary sub-coil and a second shield coil the electrical current and the second electrical current cooperating to produce a gradient coil field of view that is different from a gradient coil field of view produced without the selective removing and circulating of second electrical current.
 27. A method for designing a gradient coil for a magnetic resonance imaging apparatus, the method including: computing current densities on coaxial inner and outer cylindrical surfaces using constraints on the stored energy and gradient field linearity and uniformity, the computed current densities being generally non-zero at non-coplanar boundary contours that define coil edges; and selecting primary coil turns, shield coil turns, and coil jumps to approximate the computed current densities on the inner and outer cylindrical surfaces.
 28. The method as set forth in claim 27, wherein the computing of current densities includes: further constraining the computing to produce a substantially spatially constant eddy current effect produced by the eddy currents in a cylindrical cold shield that surrounds the outer cylindrical surface. 